Evolution of two-gap behavior of the superconductor FeSei_ 
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The superfluid density, p s , of the iron chalcogenide superconductor, FeSei-^, was studied as a 
function of pressure by means of muon-spin rotation. The zero-temperature value of p 3 increases 
with increasing transition temperature T c (increasing pressure) following the tendency observed for 
various Fe-based and cuprate superconductors. The analysis of p 3 (T) within the two-gap scheme 
reveals that the effect on both, T c and p s (0), is entirely determined by the band(s) where the large 
superconducting gap develops, while the band(s) with the small gap become practically unaffected. 
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Since the discovery of Fe-based high-temperature su- 
perconductors (HTS) much effort is devoted to the in- 
vestigation of their superconducting mechanism. While 
some properties of Fe-based HTS are reminiscent of the 
cuprate HTS (as, e.g., their layered structure, the prox- 
imity to a magnetic phase, the universal "Uemura" scal- 
ing between the superfluid density, p s , and the transi- 
tion temperature, T c ), the differences between both com- 
pounds families are much more remarkable. Hence, the 
superconductivity in Fe-based HTS originates within the 
d-orbitals of the Fe ion, which are normally expected to 
lead to pair-breaking effects [lj. For the Fe-based HTS, 
several disconnected Fermi-surface sheets contribute to 
superconductivity, as revealed by angle-resolved photoc- 
mission spectroscopy 0,0] ■ Furthermore, indications for 
multi-gap superconductivity was obtained from tunnel- 
ing magnetic torque 0], point contact [?} and in- 
frared spectroscopy [|| experiments, as well as from spe- 
cific heat [§] , first and second critical field OH [HI , and 
superfluid density 
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studies. The multi- 
gap superconducting state positions the Fe-based HTS 
together with MgB2 - the most famous double-gap su- 
perconductor discovered to date. However, it is worth 
mentionning that in spite of the fact that the two-gap su- 
perconductivity was detected for Fe-based HTS belong- 
ing to different families (as e.g. 1111: [3, @, [ll|, 0j; 
122: @, i, 0, & M, El El; Oil: 0) a systematic 
study of this phenomenon within one given family was 
not yet performed. In this paper we report on the evo- 
lution of two-gap behavior in the iron chalcogenide su- 
perconductor FeSei-a;. The transition temperature was 
changed within the range 8.3 < T c < 12.8 K by apply- 
ing an external pressure p between and 0.84 GPa. At 
each particular pressure the superfluid density p s was 
obtained from the in-plane magnetic penetration depth 
A~ b 2 (T) oc p s studied by means of muon-spin rotation, 
^SR. The analysis of A~ fc 2 (T,p), performed by solving self 
consistently the gap equations derived within the two- 
gap scheme 0, [I3]j reveals that the main effect on T c (p) 



and A~ fc 2 (T, p) tx p s (T,p) arises from the energy band(s) 
where the large superconducting gap, Ai, develops. The 
zero-temperature values of Ai, the contribution of this 
gap to the superfluid density \~ bl , as well as the effec- 
tive coupling constant An increase almost linearly with 
increasing T c (increasing pressure). In contrast, the con- 
tribution of the small gap and thus A2, A~ 6 2 2 , and A22, 
is practically pressure independent. Our results imply, 
therefore, that the transition temperature in FeSci^a; is 
entirely determined by the intraband interaction within 
the band(s) where the dominant gap is opened. 

The sample with the nominal composition FeSeo.94 was 
prepared by solid state reaction similar to that described 
in Refs. EH, [H, 0. Powders of minimum purity 99.99% 
were mixed in appropriate ratios, pressed and sealed in 
a double-walled quartz ampoule. The sample was heated 
up to 700°C followed by annealing at 400°C [IJ. The 
pressure was generated in a CuBe piston-cylinder type 
of cell especially designed to perform piSK experiments 
under pressure [2l| . As a pressure transmitting medium 
7373 Daphne oil was used. The pressure was measured 
in situ by monitoring the pressure shift of the supercon- 
ducting transition temperature of Pb and/or In. The 
/iSR experiments were carried out at the pEl beam line, 
Paul Scherrer Institute, Switzerland. Zero-field (ZF) and 
transverse-field (TF) /iSR experiments were performed at 
temperatures ranging from 0.24 to 50 K. For TF measure- 
ments the external magnetic field p>oH = 10 mT was ap- 
plied perpendicular to the muon-spin polarization. Typ- 
ical counting statistics was ~ 5 — 7 • 10 6 positron events 
for each data point. 

The results of the ZFuSR experiments were previ- 
ously reported in Rcf. 22. It was shown that up to 
p ~ 0.8 GPa the ZF response of FeSeo.94 is determined 
by the contribution of the dilute Fe moments, in analogy 
with what was observed for ambient pressure measure- 
ments of FeSeo.ss [HI- At p = 0.84 GPa static mag- 
netism was found to occupy approximately 10% of the 
sample volume at T ~ T c and it decreases down to ~ 5% 
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FIG. 1: (Color online) TF pSR time-spectra (fi H = 10 mT) 
of FeSeo.94 measured below (T — 0.24 K) and above (T = 
20 K) the superconducting transition temperature (T c ~ 
13 K) at p = 0.76 GPa. The stronger damping in the super- 
conducting state is due to the formation of the vortex lattice. 



atT- 0.25 K [22|. 

Figure [1] shows the TF /iSR time-spectra measured at 
p = 0.76 GPa above (T = 20 K) and below (T = 0.24 K) 
the superconducting transition temperature (T c ~ 13 K). 
The stronger relaxation of the muon-spin polarization at 
0.24 K relative to 20 K is due to the formation of the 
vortex lattice at T < T c . The TF /iSR data were analyzed 
by using the functional form: 



A(t) = As{t)+A PS (t) 



A s ,o e- At e-^'l' 1 cos( 7ai B s t 



-At p -<y 2 s t i /2 i 

+A PS ,o e~ als * 2/2 cos( 7 „Sp S t + 0). (1) 



Here the indexes S and PS denote the sample and the 
pressure cell, respectively, Aq is the initial asymmetry, A 
is the exponential relaxation rate caused by the presence 
of diluted Fe moments [H, 7/, = 2tt- 135.5 MHz/T is the 
muon gyromagnetic ratio, B is the internal field, and <f> is 
the initial phase of the muon-spin ensemble. The Gaus- 
sian relaxation rate, aps, reflects the depolarization due 
to the nuclear magnetism of the pressure cell, while as 
represents the depolarization in the sample arising from 
the nuclear moments and from the vortex lattice (see be- 
low). Each set of TF ^SR data taken at constant pressure 
was fitted simultaneously with As,o, Aps.o, Bps, aps, 
A, and <j>, as common parameters, and Bs and as as 
individual parameters for each temperature point. The 
exponential relaxation rate A was assumed to be temper- 
ature independent in accordance with the results of ZF 
/iSR experiments [22| . 

In an anisotropic powder sample the magnetic penetra- 
tion depth A can be extracted from the Gaussian relax- 
ation rate a sc (T) = [a 2 s {T) - al m ] 1/2 oc 1/A 2 (T), which 
probes the second moment of the magnetic field distribu- 
tion in a superconductor in the mixed state 0, 0, 24 1 . 
Here a nm is the nuclear moment contribution measured 
at T > T c . a sc can be converted into A Q {, via 
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FIG. 2: (Color online) Temperature dependence of \~ b cx p s 
of FeSeo.94 measured at p = 0.0, 0.28, 0.42, 0.58, 0.76, and 
0.84 GPa. The solid and the dashed lines are the theoretical 
curves obtained within the framework of the two-gap model 
described in the text. The insets show the temperature de- 
pendences of the large (Ai) and the small (A2) gap. 



where $0 = 2.068 • 10 15 Wb is the magnetic flux quan- 
tum. The measured A~ b 2 (T,p) of FeSeo.94 at p = 0.0, 
0.28, 0.42, 0.58, 0.76, and 0.84 GPa are shown in Fig. [2] 

The experimental A~ b 2 (T) data were analyzed by us 
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ing the two- gap model presented recently in Rcfs 
Following [16j , the coupled gap equations for a supercon- 
ductor with fc-dependent energy gaps: Ai, A2, intraband 
pairing potentials: Vn, V22, and interband interaction 
potentials: V12, V21 axe determined as: 



x tanh ■ 



2k B T 



(3) 



Wu = 0-00126 $ 2 /A Q 4 b , 



(2) 



Here i = 1, 2 is the band index and the sums are taken 
within the corresponding energy bands in fc-space. Con- 
sidering isotropic s-wave gaps [Aj(fcj) = Ai] and ne- 
glecting the momentum dependence of the band energies, 
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Eq. ([3]) reduces to: 
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(4) 

Here Lo Bi is the phonon cutoff frequency (Debye fre- 
quency; note that these cutoffs are expected to be dif- 
ferent for both bands u>d 1 ^d 2 an< ^ -A^(O) 1S the 
partial density of states at the Fermi level. For conve- 
nience the sums were also converted into integrals 
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A further simplification of Eq. (J4j) can be made by us- 
ing the notation of the coupling constant Ay = A;(0)Vy 
introduced by Kogan et al. 25[ and assuming similar 
cutoff frequencies for both bands (uJdi = ^>d 2 = ^d)- 



E 
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E 2 + A 2 



E 2 + A 2 
: tanh 1 — -dE. 



2k B T 



(5) 

The advantages to use the above equation rather than 
Eq. ((4)) is that (i) within the notation of Kogan et al. 
[251 ] A12 = A21 and (ii) the total number of parameters 
needed to evaluate Ai(T) and A 2 (T) reduces from 8 in 
case of Eq. (01 to 4 in case of Eq. ([5]). 

With the known Aj(T) and A 2 (T), 
taincd by decomposing it into two components A~^, and 



A . can be ob- 
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\ ab 2 so that: 
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Here / = [1 + cxp(E/kBT)]~ 1 is the Fermi function. 

The analysis of A~ b 2 (T, p) by using the above described 
model was made by solely evaluating A22, A~ 61 (0), and 
A~ b 2 2 (0). The parameters An, A 12 , and lod were taken 
as follows: 

Ai 2 : Our numerical analysis reveals that the step- like 
change of A~ b 2 (T) at T ~ 2.5 K (see Fig. requires 
the interband coupling constant A12 to be very small 
(A12 ~ 10~ 3 or smaller). This implies that the band(s), 
where the large and the small superconducting energy 
gaps are open, become only weakly coupled. Note that a 
similarly small interband coupling constant was obtained 
by Kogan et al. for the superconductor VsSi 25 1. 
An: The fact that the interband coupling in FeSei_ K is 
weak thus suggests that the transition temperature T c is 
mainly determined by the coupling within the band(s) 
where the large superconducting gap is opened. As- 
suming that the larger gap is Ai, T c is defined when 
Ai(T) = 0, so that according to Eq. ©: 
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FIG. 3: (Color online) (a) Dependence of \~£(T = 0) on 
T c . The inset is the "Uemura relation" for Fe-based HTS 
with some of data obtained to date (see Ref. [Til . HH I23I I30I ). 
(b), (c), and (d) Dependence of the parameters obtained in 
analysis of \~^(T,p) within the framework of two-gap model 
(see text for details). 



With T c (p) measured independently (see Ref. [22j) 
Eq. (HI) allows one to obtain the value of the intraband 
coupling constant An for each particular pressure. 
lu b : The ambient pressure value of the cutoff phonon 
frequency (Debye frequency) ojd{p — 0) ~ 40 meV was 
taken from Ref. [27j . The increase of u>d with increasing 
pressure was assumed to follow: 



wd(p) =wd(0)(1 + 7p/B), 



(8) 



which is the consequence of the Gruneisen equation 
7 = — dlnujD I dhiV (7 is the Griinciscn parameter, 
B ~ 31 GPa is the bulk modulus [3], and V is the sample 
volume). The Gruneisen par ameter was assumed to be 
7 « 1 in analogy with Ref. |2_Sj. We should also emphasize 
that the parameters of the above described model are not 
very sensitive to the exact value of lo b . As an example, 
the increase (decrease) of ojd by a factor of 2 leads to 
a corresponding decrease(increase) of An obtained from 
Eq. ([5]) by ~ 15(18)%. This make our assumption about 
using a similar cutoff phonon frequency for both bands 
[udj = uj£> 2 , see Eq. (JSJ)] to be rather reliable. Note that 
similar conclusion was also reached by Kogan et al. [2^ ] . 

The parameters obtained from the analysis of A~ b 2 (T) 
by means of the model described above are summarized 
in Table |TJ The value of the interband coupling constant 
A12 = 0.005 was kept fixed. The red and the blue dashed 
lines in Fig. [2] correspond to the contribution of the large, 
A~ 61 , and the small, A~ b 2 2 , superconducting gaps to the 
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TABLE I: Summary of the pressure studies of FeSei-z. The meaning of the parameters is: p - pressure; T c - transition 
temperature; cup - Debye frequency; A12 - interband coupling constant, A11/A22, Ai(0)/A2(0), A~, 2 \(0) / 'A~ 6 S ! 2 (0) - intraband 
coupling constant, zero-temperature value of the gap, zero temperature value of superfluid density component within the band 
1/2, respectively. 

p T c cod An A22 A 12 Ai(0) A 2 (0) A^(0) A^Of 

(GPa) (K) (meV) (meV) (meV) (fim~ 2 ) (/im~ 2 ) 

0.00 8.3(1) 40.0 0.241(1) 0.195(1) 0.0005 1.27(1) 0.487(6) 3.67(25) 3.53(25) 
0.28 10.2(1) 40.36 0.253(1) 0.195(1) 0.0005 1.56(1) 0.494(6) 5.19(23) 3.61(23) 
0.42 11.6(1) 40.54 0.261(1) 0.195(1) 0.0005 1.77(1) 0.498(6) 5.20(20) 4.00(18) 
0.58 12.4(1) 40.74 0.265(1) 0.194(1) 0.0005 1.89(1) 0.491(6) 5.45(19) 4.11(20) 
0.76 12.8(1) 40.98 0.267(1) 0.195(1) 0.0005 1.94(1) 0.504(6) 6.44(29) 4.12(25) 
0.84 12.6(1) 41.08 0.266(1) 0.193(1) 0.0005 1.92(1) 0.480(6) 6.74(29) 3.47(29) 



total superfluid density, solid lines. The temperature de- 
pendences of the large, Ai, and the small, A2, gaps are 
shown in the corresponding insets. 

In order to check how the change of T c affects the 
energy bands where the large and the small supercon- 
ducting gaps are supposed to be open, we plot in Fig. [3] 
the parameters A" 2 (0), A^O), A" 2 2 (0), Ai(0), A 2 (0), 
An = iVi(0) Vn, and A 22 = N 2 (0) V21 as a function 
of the transition temperature T c . From the obtained 
data the following conclusions can be drawn: (i) The 
zero-temperature superfluid density p s (0) oc A~ b 2 (0) in- 
creases with increasing T c thus following the "Uemura" 
relation established recently for various Fe-bascd HTS, 
see Fig. |3] (a) and Refs. 0, OH 0, |H (ii) The electronic 
bands, where the large and the small gap are opened, are 
affected by the pressure quite differently. The increase 
of T c with pressure leads to an almost linear increase of 
the superfluid density component A~ b 2 1 (0), the supercon- 
ducting energy gap A 1,0 as well as the effective coupling 
constant An = iVi(0) V\\. On the other hand, both A 2; o 
and A 22 = A^ 2 (0) V22 stay almost constant, while A~ b 2 2 (0) 
increases with increasing T c only slightly [see Figs. [3j (b), 
(c) and (e)]. Bearing in mind that the "large gap" and 
the "small gap" bands are only weakly coupled (the in- 
terband coupling constant Ai 2 is estimated to be of the 
order of 5 • 10~ 4 or less, see Table Q} one may conclude 
that in the range of < p < 0.84 GPa the pressure effect 
on both T c and A~ b 2 is solely determined by the bands 
exhibiting the large superconducting gap. 

To conclude, the superfluid density p s oc A~ b 2 was stud- 
ied as a function of pressure and temperature in the su- 
perconductor FeSei_ 2 - by means of /iSR. The analysis of 
p s (T) within a two-gap scheme reveals that the effect on 
both, T c and p s (0), is entirely determined by the band(s) 
where the large superconducting gap develops. Our re- 
sults suggests that for 011 family of Fe-bascd HTS the in- 
traband interaction is most probably the leading pairing 
interaction determining the superconducting properties. 
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